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Mode transformation and frequency change with height in 3D numeri- 
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Three-dimensional numerical simulations of magnetoacoustic wave propagation are performed in a sunspot atmosphere 
with a computational domain covering from the photosphere to the chromosphere. The wave source, with properties 
resembling the solar spectrum, is located at different distances from the axis of the sunspot for each simulation. These 
results are compared with the theory of mode transformation and also with observational features. Simulations show that 
the dominant oscillation frequency in the chromosphere decreases with the radial distance from the sunspot axis. The 
energy flux of the different wave modes involved, including de Alfven mode, is evaluated and discussed. 
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1 Introduction 

Observations of waves in sunspots show a variety of be- 
haviours depending on the height and the region observed. 
In the umbral photosphere their power spectrum has a peak 
at 3.3 mHz (period of 5 min) and their amplitudes are a- 
round a hundred m s~^. The frequency of the power peak 
varies with the height, and at the chromosphere the oscilla- 
tory pattern is dominated by waves in the 5-6 mHz band (pe- 
riod of 3 min) with amplitudes of several kilometers per sec- 
ond which present a saw-tooth profile. Several mechanisms 
have been proposed to explain the chromospheric three min 



Most of the simulations of mode refraction, transfor- 
mation, etc. have been developed for waves with high fre- 
quencies, i.e. above the cut-off frequency of the solar atmo- 
sphere, and using 2D atmospheres. We aim to study these 
phenomena through numerical simulations of waves with 
the periods observed and using realistic atmospheres of sun- 
spots in a 3D case. One of the most important benefits of 
3D simulations of wave propagation in magnetized atmo- 
spheres is the app earance of the Alfven mode. The work by 



Felipe et al. I (1201 0) points out that the conversion from the 



slow to the Alfven mode is effective when the wave vec- 
tor forms a certain angle different from zero with the mag- 
ute os cillations: a resonant chromospheric cavity ( Zhugzhdaet afl^.f ic field. C onversion to Alfven waves only was produced 
1985 | ) : non-linear interaction of photospheric modes ( Gurman^iiEM^er was located out of the axis of the sunspot, 

where the magnetic field has some inclination. In this con- 
tribution, we explore further the conclusions of that work 
to study waves with 3-5 min periodicity found in sunspots, 
and analyze several simulations with the driver located at 
different distances from the center of the sunspot. 



19841) : and slow magneto-acoustic mode field aligned prop- 



agation fro m the photosphere to the chromosphere in the 5-6 
mHz band dCenteno et al.ll2006l) . 



The interpretation of the observed oscillations in terms 
of MHD waves possess several difficulties. The velocity of 
propagation of the magneto-acoustic modes depends on the 
acoustic velocity (c^) and Alfven velocity (va), both strat- 
ified with height. The magnetized atmosphere of sunspots 
changes from being gas dominated in the photosphere and 
below, (c5 > va) to field dominated at the high chromo- 
sphere {cs < va)- Around the layer where cs = va the 
phase speeds of all modes are similar and different waves 



2 Numerical simulations 



We have developed a 3D non-linear MHD code (IFelipe et al. [ 
2010) to study the wave interaction with magnetic struc- 

can interact, whic h produces the mode tran sformation (ij^^n p^^fl ^|^^ excite a MHS model of sunspot, constructed fol- 

I2003l :[c^l200ai]aome^^^ & Colladosll2006l) . The direc- ^^^mg Kh omenko & Collados| (|2008|), with a perturbation 

tion and efficiency of the mode transformation depends on 

the frequency of the wave an d the angle b etween the wave 

vector and the magnetic field (ICallyll2006l) . 
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which drives a wave spectrum close to the solar one. The 
computational domain covers 60 Mm in both horizontal di- 
mensions, with the axis of the sunspot located at the center, 
and spanning from z = —2.5 to z = 1 Mm. The height 
z = Mm corresponds to the layer where t^qq = 1 at quiet 
Sun photosphere. The full set of simulations consists of 3 
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Fig. 1 Variations of longitudinal velocity scaled by a factor p^l^ at an elapsed time t=820 s after the beginning of the 
simulations with the driver force at 5 Mm from the axis of the sunspot. Grey inclined lines are magnetic field lines. 




Fig. 2 Variations of transversal velocity scaled by a factor p^/^ at an elapsed time t=820 s after the beginning of the 
simulations with the driver force at 5 Mm from the axis of the sunspot. Grey inclined lines are magnetic field lines. 



runs with the driver located at Xq = 5,10 and 15 Mm from 
the axis, respectively, and lasting 30 min of solar time. 



Figures [T] and [2l show the longitudinal and transversal 
components of the velocity, respectively, of the simulation 
with the driver located at 5 Mm from the axis of the sunspot. 
The driver mainly excites a fast acoustic mode below the 
(3 = 1 layer. When it reaches this layer it is transformed 
into a slow acoustic mode (visible in blue in high layers at 
X = — 5 Mm in the longitudinal velocity. Fig. [T]), which 
propagates upwards to the chromosphere; and a fast mag- 
netic mode (visible in the transversal velocity, Fig.O, which 
is reflected back to the photosphere due to the gradients of 
the Alfven speed. When it reaches the /3 = 1 layer, it is 
transformed again into a fast acoustic mode and a slow mag- 
netic mode in the high-/3 region, which appears in Figs. [T] 
and O respectively, as concentric rings around and below 
z = Mm. There is also conversion to the Alfven mode, 
although it is not visible in these figures. We will discuss it 
in Section m 



3 Frequency change with height 

Figure [3] shows the power spectra of the longitudinal ve- 
locity at two different heights at the location of the driver 
for the simulation with Xq = 10 Mm. The blue line corre- 
sponds to the photospheric height z = —1 Mm, where the 
driver was imposed, while the red line is the power spectra 
at z = 1 Mm, i.e., at the chromosphere. The photospheric 
power spectra peaks at 3.3 mHz (5 min) and it has a sec- 
ondary peak around 5.4 mHz. The maximum contribution 
to the power at the chromosphere is at 5.3 mHz (3 min), 
so the simulation reflects the observed change of frequency 
with height. 

The power spectra at the chromosphere in Fig. [3] shows 
some secondary peaks, which correspond to the harmon- 
ics of the 3 min signal, due to the importance of the non- 
linearities at these high layers, where the longitudinal ve- 
locity reaches an amplitude of 4.5 km s~^. Waves with fre- 
quecies below the cut-off (dashed vertical line) form evanes- 
cent waves, which can not propagate energy upwards, while 
waves with frequencies above the cut-off (as the secondary 
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Fig. 3 Normalized power spectra of the longitudinal ve- 
locity at the photosphere (blue line) and the chromosphere 
(red line) at the location of the driver for the simulation with 
Xq = 10 Mm. Asterisks mark the measured values. The 
vertical dashed line marks the maximum value of the cut- 
off frequency at this position. 
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Fig. 4 Height variation of the dominant frequency of os- 
cillations for the 3 simulations with the driver located at 
5 (red line), 10 (green line) and 15 (blue line) Mm of the 
sunspot axis. 



peak of the photospheric power spectra) do propagate en- 
ergy upwards. Their amplitude increases according to the 
drop of the density. 

Figure [H shows the frequency of the dominant power 
peak at all atmospheric heights. The frequency of the peak at 
the chromosphere decreases with the radial distance, since 
the cut-off frequency is higher near the center of the sun- 
spot. At photospheric deeper layers the dominant frequency 
is around 3.3 mHz (5 min) for all simulations with different 
locations of the driver, increasing abruptly around z = 
Mm (depending on the case), when the power at frequencies 
above the cut-off becomes more important than the evanes- 
cent 3 mHz oscillations. 
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Fig. 5 Acoustic (top) and magnetic (bottom) flux for the 
simulation with Xq = 5 Mm averaged over the station- 
ary stage of the simulations in logaritmic scale. Horizontal 
white line is the height where sound velocity and Alfven 
velocity are equal. 
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Fig. 6 Acoustic (top) and magnetic (bottom) flux for the 
simulation with Xq = 10 Mm averaged over the station- 
ary stage of the simulations in logaritmic scale. Horizontal 
white line is the height where sound velocity and Alfven 
velocity are equal. 

4 Energy balance 

The driver that we introduce as a perturbation mainly ex- 
cites a fast acoustic wave in the region where sound speed 
is higher than Alfven speed. The wave energy flux can be 
calculated as 

F = Fac + F„,ag = Piv + Bi X (v X Bo)//io, (1) 

where pi and Bi are the perturbations in the pressure and 
magnetic field, respectively, v is the velocity, Bq is the equi- 
libirum magnetic field, and /io is the magnetic permeability. 
FiguresOm and[7]present the acoustic (top panels) and 
magnetic (bottom panels) energy fluxes obtained from Eq. 
([B. They clearly show that most of the energy introduced by 
the driver keeps below the layer cs = va (white horizontal 
line), since the low frequency (3.3 mHz band) slow acous- 
tic mode can only propagate horizontally, and the fast mag- 
netic mode above this layer returns towards the photosphere 
and generates new acoustic and magnetic flux after its trans- 
formation. Only the slow acoustic mode with frequencies 
above the cut-off frequency that propagates upwards along 
field lines, visible in top panel of Figs. [5] — [7] at a radial dis- 
tance around the location of the driver and z between 0.5 
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Fig. 7 Acoustic (top) and magnetic (bottom) flux for the 
simulation with Xq = 15 Mm averaged over the station- 
ary stage of the simulations in logaritmic scale. Horizontal 
white line is the height where sound velocity and Alfven 
velocity are equal. 

and 1 Mm, supplies energy to the chromosphere, since the 
energy of the Alfven mode is negligible, as we will show 
later. Most of the magnetic flux above the layer cs = va 
corresponds to the fast magnetic waves which are being re- 
flected. 

The Alfven mode is visible in Fig.[8l It shows the longi- 
tudinal magnetic energy flux calculated from Eq. ([T]) when 
V and Bi are proj ections along the polariz ation direction 
of the Alfven wave ( iCally & Goossenj 2008 ). for the driver 
located at 5 Mm from the axis of the sunspot. The top panel 
is a cut in the plane x — z aty = —1.3 Mm, and it is nor- 
malized to its maximum at every height. The bottom panel 
is a horizontal cut in the plane x — y at z = 0.5 Mm. The 
Alfven direction projections from lCally & GoossensI (|2008|) 
are valid only in asymptotic case /3 <C 1. It means that in 
the top panel we only can indentify the Alfven mode energy 
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Fig. 8 Longitudinal magnetic flux of the Alfven mode at 
t = 330 s. Top: Vertical cut in the plane y = —1.3 Mm, 
normalized at every height. Vertical white lines are mag- 
netic field lines and horizontal white line is the layer where 
flux in the layers above the fine c^ = v\, and the flux below c\ = v\. Bottom: Horizontal cut in the plane z = 0.5 Mm. 



this height is not meaningfull. The Alfven wave propagates 
upwards along field lines. The transformation from the fast 
acoustic mode to the Alfven wave is null at the plane ^ = 
Mm, where the driver is located, but out of this plane there 
is energy flux in the Alfven mode, as can be seen in the 
bottom panel. However, for the inclination of the magnetic 
field and the direction of propagation of the wave in these 
simulations, the trasformation to the Alfven mode is very 
ineffective, and the magnetic flux of this wave at the height 
z = 0.5 Mm is 40 times lower than the acoustic flux of the 
slow mode. 

5 Conclusions 

- The mechanism that produces the change in frequency 
of oscillations in the umbra from the photosphere to the 
chromosphere is the linear propagation of waves with 
3 minute power which come directly from the photo- 
sphere. 

- The dominant wave frequency in the chromosphere de- 
creases with the radial distance, due to the reduction of 
the cut-off frequency far from the axis. 



Thin dashed lines are contours of equal magnetic field, with 
its value indicated in Gauss. In both panels, black dashed 
lines mark the location of the other plot. 



- Only waves in the 5-6 mHz frequency band can pro- 
vide energy to the chromosphere through slow acoustic 
mode, since waves with lower frequencies form evanes- 
cent waves and do not carry energy, and the fast mag- 
netic mode is reflected back towards the photosphere. 

- The conversion to the Alfven mode is effective where 
the magnetic field has some inclination, but its energy 
flux is small compared to the acoustic energy flux. 

- Most of the energy is confined to the region below the 
height where the sound and Alfven velocity are similar. 
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